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ABSTRACT 


Systems  analysis  techniques  are  applied  to  the  problem  of  pilot- 
induced  aircraft  oscillations  (PIO).  Mathematical  models  used  for  the 
pilot  prior  to  and  during  PIO,  as  well  as  the  use  of  various  pilot  and 
subsystem  describing  functions,  are  reviewed.  Several  examples  of  PIO 
causes,  analyses,  and  cures  are  given,  including  linear  and  nonlinear 
longitudinal  PIO.  The  closed-loop  describing  function  for  a  rate- 
limited  position  servo  is  derived  in  the  Appendix. 


rORWORD 


This  report  summarizes  several  years  of  intermittent 
work  conducted  at  Systems  Technology,  Inc.,  on  the  prob¬ 
lem  of  pilot-induced  oscillations.  The  major  portions 
of  the  efforts  reported  were  performed  either  as  part  of 
other  handling  qualities  studies  under  Contract 

33(61 6)  -5661  for  the  Flight  Control  Laboratory  of  the 
Air  Force  Systems  Command,  or  as  part  of  an  investigation 
of  T-38A  and  general  PIO  problems  for  Northrop-Norair 
under  Lietter-Agreetaent-  8944-62-  582~dated  24  July  1 962, 
and  subsequent  amendments.  This  report  itself  has  been 
almost  totally  drafted  under  the  Norair  agreement. 

The  major  contributors  to  these  efforts,  and  to  this 
report,  are  listed  as  authors.  Other  STI  engineers  who 
have  also  contributed  in  some  way  or  other  to  STI's  PIO 
studies  include:  D.  Graham,  W.  A.  Johnson,  V.  J.  Kovacevich, 
C.  P.  Shortwell,  J.  Wolkovitch,  and  C.  D.  Wezeman.  The 
authors  are  most  grateful  to  both  Mr.  D.  L.  Hirsch  of  the 
Vehicle  Dynamics  and  Control  Branch,  Norair,  and 
Mr.  R.  J.  Wasicko,  who  served  as  the  AF  Project  Engineer — 
their  guidance,  constructive  comments,  support,  and  patience 
have  played  an  extremely  important  role  in  these  studies. 
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A  pilot-induced  oscillation  (PIO)  is  an  inadvertent  sustained 
oscillation  of  the  pilot-vehicle  system.  Throughout  aviation's  history 
PIO's  of  various  degrees  of  severity  have  appeared  sporadically.  Indeed, 
the  static  instability  of  the  Wright  Flyer  resulted  in  a  mild  longitudinal 
oscillation  of  that  most  early  pilot-vehicle  system.  With  very  high  per¬ 
formance  aircraft  having  fully  powered  control  systems,  PIO's  have  become 
more  frequent;  .and  also  more  lethal,  because  relatively  small  pilot  forces 
can  cause  a  rapid  buildup  to  catastrophic  loads.  Many  of  /these  cases 
appear',  embarrassingly,  during  final  flight  tests  or  early  production 
stages;  thus  there  are  few  well-documented  cases  on  public  record.  Some 
exceptions  are  reported  in  Refs.  1-4,  18,  and  19*  Two  of  these  references 
(l  and  4)  present  a 'rare,  and  unusually  complete,  time  history  of  a  PIO. 
This  dramatic  recording,  which  has  great  value  in  providing  physical 
appreciation  for  the  phenomena  involved  in  PIO's,  is  reproduced  in  Fig.  1. 
The  PIO  occurred  in  an  early  verson  of  the  T-58  aircraft  (since  modified 
to  completely  eliminate  the  possibility  of  recurrence)  following  the  shut¬ 
off  of  an  oscillating  pitch  damper  in  a  mistrimmed  position.  Besides  pro¬ 
viding  a  beautiful  example  of  events  in  a  PIO  buildup  sequence  which  will 
be  elaborated  further  below,  this  record  with  its  ±8g  oscillations  offers 
graphic  experimental  evidence  of  a  situation  which  pilots  will  not 
willingly  duplicate. 

The  often  violent  and  unexpected  nature  of  many  forms  of  PIO's  has 
usually  lead  to  a  prompt  attempt  to  "cure"  the  condition  (quietly). 
Although  there  is  no  generally  valid  theory  for  predicting  PIO's,  two 
lines  of  approach  (Ref.  2)  can  be  indicated: 

1 .  Compare  the  open-loop  dynamic  characteristics  of  the  aircraft 
plus  manual  controls  with  those  dynamic  characteristics  of 
other  craft  which,  as  shown  by  experience,  are  compatible 
with  the  presence  or  absence  of  PIO's. 
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2.  Perform  closed- loop  stability  analyses  of  the  pilot-vehicle 
system  using  appropriate  describing  function  representations 
for  the  pilot,  vehicle,  and  manual  control  system. 

At  one  limit  the  first  approach  leads  to  a  catalog  of  vehicle-plus- 
control-  system  sinusoidal-input  describing  functions  which  represent  either 
PlO-prone  or  PlO-resistant  configurations.  By  their  nature  such  data  are 
extremely  valuable  as  points  of  departure,  although  they  are  usually  too 
specific  to  be  directly  useful  on  other  vehicle  configurations.  An  attempt 
to  generalize  on  PIO-  resistant  characteristics  for  a  linear  vehicle  is 
given  in  Ref.  2.  Referring  to  the  airframe ' 8' frequency  response  for 
dynamic  stick  force  per  g,  it  is  noted  that:  [when  the  steady-state  stick 
force  per  g  is  itself  satisfactory]  . .experience  shows  that  unsatisfac¬ 
tory  dynamic  characteristics  of  the  pilot-airframe  combination  and  pilot- 
induced  oscillations  are  avoided  when  the  amplitudes  of  the  control-force- 
gradient  frequency  response  are  never  less  than  the  low-frequency  asymptote." 
In  practice  this  corresponds  to  a  short-period  damping  ratio  of  1  tyf*  or 
higher.  More  recent  experience  (e.g.,  Refs.  5  and  6)  indicates  that  such 
a  high  danping  ratio  is  not  essential  to  avoid  PIO’s,  so  the  quoted  state¬ 
ment  represents  a  sufficient  rather  than  necessary  condition.  However, 
some  such  specification  is  desirable  for  its  simplicity,  and,  perhaps,  as 
still  more  experience  is  gained,  a  similar  statement  can  ultimately  be 
refined  to  include  both  necessary  and  sufficient  conditions. 

In  the  meantime  one  can  either  overspecify  on  the  basis  of  applicable 
experience  (the  first-listed  approach)  or  treat  PIO  problems  using  a 
servoanalytic  attack  (the  second- listed  approach).  The  latter  approach 
ha6  been  used  for  some  time  with  considerable  success  to  provide  physical 
explanations  of  PIO  occurrences  and  to  enhance  appreciation  of  the  effects 
of  fixes  on  PIO  tendencies  (e.g.,  Refs.  1-4,  '(  and  20).  The  techniques  used 
provide  a  measure  of  prediction  as  well  as  understanding.  The  approach  is 
limited  by  incomplete  knowledge  of  the  range  and  applicability  of  possible 
human-pilot  describing  functions.  However,  the  state  of  such  knowledge 
has  progressed  through  the  years,  especially  with  respect  to  closed-loop 
handling  qualities  theory,  so  some  statements  can  now  be  made  with  greater 
certainty  than  previously.  One  purpose  of  this  report  is  to  bring  the 
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state  of  the  art.  of  PIO  technology  up  to  date  via  the  systems  analysis 
approach  which  has  been  successful  in  the  related  field  of  handling 
qualities  research  (Refs.  9— 1 2) • 

The  over-all  system  to  be  analyzed  in  PIO  problems  encompasses  complex 
dynamic  interactions  between  the  pilot,  artificial  feel  devices,  powered 
control  surfaces,  airframe,  and  display  subsystems.  To  this  general  com- 
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plexity  are  added  the  following  complications: 

Adaptive  behavior  by  the  pilot 

The  presence  of  unavoidable  or  intentional  nonlinearities 
(such  as  friction,  dead  zones,  and  preload)  in  the  control 
system 

« 

Often  neglected  aerodynamic  terms  (such  as  the  elevator's 
contribution  to  the  normal  acceleration  felt  by  a  pilot 
and,  for  lateral  PIO,  the  aileron's  adverse  yaw) 

In  turn,  describing  the  adaptive  pilot  requires  a  valid  mathematical  model, 
the  nonlinearities  require  consideration  of  the  specific  task  and  disturb¬ 
ance  inputs,  and  the  aerodynamic  subtleties  require  a  thorough  understand¬ 
ing  of  the  airframe  transfer  functions.  In  spite  of  the  total  complexity 
of  the  pilot-vehicle  system  containing  these  elements,  PIO's  can  and  have 
been  successfully  analyzed  and  predicted  by  careful  system  analyses  of 
the  types  tc  be  described  herein. 

The  report  starts  with  a  brief  resume  of  the  systems  approach  to  the 
analysis  of  the  pilot-vehicle  system  in  Section  II.  This  lays  the  founda¬ 
tion  for  Section  III,  which  introduces  and  categorizes  the  basic  causes  of 
PIO.  Section  IV  contains  selected  examples  of  the  system  analysis  of  PIO 
situations,  and  covers  both  the  diagnostic  and  the  curative  aspects  of  the 
problem  encountered.  Mathematical  derivation  of  a  describing  function 
needed  in  Section  IV,  but  not  available  elsewhere,  is  contained  in  the 
Appendix. 


4 


MOTION  II 

SHE  PILOT-VEHICLE  SYSTEM 

A,  OENSBAL  DESCRIPTION 

Among  the  significant  attributes  of  a  human  pilot  are  his  ability  to 
establish  a  wide  variety  of  pilot-vehicle  system  organizations  (i.e., 
many  different  system  structures  or  effective  block  diagrams)  and  adjust¬ 
ments  therein.  This  repertory  of  behavior  is  so  extensive  that  the  pilot, 
as  an  adaptive  controller,  has  capabilities  which  exceed  those  of  the  most 
sophisticated  unmanned  control  system.  From  a*  systems  analysis  standpoint 
this  variety  may  at  first  seem  discouraging.  For  many  flight  control 
situations,  however,  further  consideration  indicates  that  all  is  not  hope¬ 
less.  In  controlling  any  complex  system,  successful  behavior  is  very 
narrowly  limited.  The  very  nature  of  the  requirements  for  "good"  control 
system  performance  and  the  restrictions  imposed  by  the  dynamic  character¬ 
istics  of  the  controlled  element  constrain  the  successful  human  pilot  to 
operate  in  accord  with  will-established  "laws."  When  well  trained  and 
motivated,  or  when  the  imposed  restrictions  are  severe,  the  performance 
of  the  pilot  and  the  system  can  be  predicted  with  an  accuracy  sufficient 
for  many  engineering  purposes.  The  prediction  is  both  qualitative  and 
quantitative.  In  short,  the  man-machine  system  can  be  made  amenable  to 
mathematical  analysis  (Refs.  1  >,  and  16). 

Experiment  has  shown  that  the  system  organizations  and  adjustments 
adopted  by  the  pilot  for  single-loop  or  uncoupled  multiloop  feedback 
systems  are  consistent  vrith  those  of  "good"  feedback  systems  in  general. 
This  presumes,  of  course,  that  the  adjustments  which  are  required  do  not 
violate  the  capabilities  of  the  hianan  operator.  Thus,  the  system  will  be 
organized  and  adjusted  as  if  done  by  a  "super-servoanalyst. "  In  the  case 
of  more  complicated  multiloop  systems,  direct  experimental  data  on  pilot 
adjustment  and  organization  do  not  yet  exist.  However,  comparison  of 
experimental  results  with  analytical  inferences  obtained  with  an  extended 


pilot  model  strongly  suggests  that  the  pilot  model  and  analysis  methods 
have  a  continuing  practical  validity.  In  fact,  as  the  number  of  possible 
display  and  control  alternatives  increases,  the  most  fruitful  criteria 
yet  found  for  choosing  the  dominant  loop  closures  are  simply  those  which 
the  "super-servoanalyst"  would  choose. 

B.  XBPUP-OOTFUE  MODEM 

The  models  used  in  pilot-vehicle  systems  analyses  are  input-output 
operators,  that  is,  transfer  functions  and  describing  functions.  In  the 
case  of  constant-coefficient  linear  systems  the  transfer  function  is 
simply  given  by  the  Cramer's  Rule  solution  of  the  equations  of  motion 
when  written  in  Laplace  transform  form.  Thus  transfer  functions  are 
ratios  of  rational  polynomials  in  the  Laplace  operator,  s;  the  denominator 
is  the  characteristic  equation  of  the  system,  while  the  numerator  connects 
a  specific  output  with  specific  forcing  terms.  Substitution  of  ja>  for  s 
in  the  transfer  function  will  yield  the  "frequency  response."  The  frequency 
response  relates,  via  an  input/output  amplitude  ratio  and  phase  angle 
difference,  tn-..  fundamental-sinusoid  component  of  the  system  output  to  the 
sinusoidal  input.  (If  the  system  is  stable  the  total  output  will  approach 
the  sinusoidal  component  of  the  output  as  time  goes  by  and  the  transients 
die  away) . 

The  extension  of  the  basic  notion  of  an  input— output  operator  to 
describe  the  input-output  characteristics  of  nonlinear  elements  is  fraught 
with  subtle  mathematical  difficulties  (Ref.  15)  but  can  easily  be  visualized. 
Physically,  in  most  oscillatory  circumstances  the  input  to  the  nonlinearity 
is  usually  close  to  a  sinusoid,  whereas  the  output  of  zhe  nonlinearity  is 
periodic  but  may  depart  considerably  from  a  sinusoidal  shape.  This  distorted 
waveform  reflects  the  presence  of  higher  harmonic  components  induced  by 
the  nonlinearity.  In  their  transmission  about  the  loop,  these  harmonics 
are  substantially  suppressed  in  amplitude  by  the  lags  of  the  airframe  and 
control  system.  Because  they  largely  disappear  in  transit,  the  input  to 
the  nonlinearity  is  scarcely  affected  by  these,  higher  harmonics,  and  thus 
retains  •'ts  nearly  sinusoidal  nature.  So,  insofar  as  stability  consider¬ 
ations  are  concerned,  the  inputs  and  outputs  of  most  of  the  nonlinear 


elements  of  Interest  in  PIO  (nonlinear  gearing,  friction,  breakout,  etc.) 
can  be  approximated  by  a  pair  of  sine  waves  of  suitable  amplitude  and 
phase  shift. 

The  appropriate  sinusoidal  approximations  for  both  the  input  and  the 
output  of  the  nonlinearity  are  the  fundamental  components  of  the  Fourier 
series  representation  of  the  actual  waveforms.  Then,  analogous  to  trans¬ 
fer  funqtions,  the  "sinusoidal- input  describing  function"  is  defined  as 
the  amplitude  ratio  and  phase  shift  between  the  fundamental  components  of 
the  input  and  output.  However,  the  output  waveform  from  a  nonlinear 
element  at  any  given  input  frequency  may  vary  as  the  input  amplitude 
varies,  so  each  nonlinearity  requires  a  family  of  describing  functions  at 
specified  input  levels.  In  fact,  it  is  the  often  adverse  changes  with 
input  amplitude  (in  the  describing  function  amplitude  ratio  and  phase 
shift)  which  lead  to  sustained  closed-loop  system  oscillations.  The  use 
of  rational  polynomials  to  represent  sirusoidal  describing  functions  for 
nonlinear  systems  is  only  justified  in  the  frequency  domain;  and  the  sub¬ 
stitution  therein  of  the  Laplace  operator,  s,  for  jcn  is  not  meaningful. 

Thus,  it  is  not  strictly  correct  to  plot  a  conventional  root  locus  diagram 
for  a  loop  closed  around  a  nonlinearity,  though  in  many  cases  this  is  done 
to  help  describe  qualitatively  the  effect  of  the  nonlinearity.  The  only 
strictly  valid  portion  of  such  nonlinear  root  loci  is  on  the  jo>  axis,  and 
cm  the  particular  locus  family  corresponding  to  a  certain  input  amplitude 
at  the  nonlinear  element.  It  also  follows  that  the  transient  response  of 
nonlinear  systems  cannot  be  calculated  from  information  in  the  sinusoidal 
describing  function. 

0.  PILOT  MODUS 

Two  phases  or  types  of  manual  control  system  behavior  must  be  considered 
in  a  PIO  analysis.  These  are: 

1 .  Prior  to  a  PIO  the  pilot  is  exerting  on -the -average 
control  of  more-or-less  random  inputs  and  has  adopted 
a  quasi-stationary  set  of  feedbacks  (displayed  or 
sensed  quantities)  and  equalisations  (gains,  lead, 
lag)  which  are  coapatible  with  "good"  control  (i.e., 
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small  error,  stability,  low  effort,  etc.)*  Several 
control  loops  may  be  present,  although  only  one  or 
two  are  usually  dominant. 

2.  After  a  PIO  is  developed,  the  apparent  airframe 
motions  change  from  a  random-like  to  a  nearly  sinu¬ 
soidal  form. 

Corresponding  to  these  two  types  of  system  behavior  are  two  quite  different 
varieties  of  pilot  and  nonlinear  element  describing  functions.  For  the 
first  phase,  at, least  some  basic  portions  of  the  total  control  system,  as 
organized  by  the  pilot,  are  compensatory  (the  total  control  system  may  have 
pursuit  elements  inserted  by  the  pilot,  see  Ref.  14,  but  these  will  not 
affect  those  system  aspects  relating  to  stability),  and  the  appropriate 
models  are  Gaussian-input  describing  functions.  For  inanimate  elements 
these  are  known  only  for  the  simplest  of  nonlinearities,  but  for  the  pilot 
the  Gaussian-input  describing  function  for  compensatory  situations  is 
fairly  well  established.  Figure  2  contains  a  brief  summary  and  Ref.  16 
a  simple,  more  complete  treatment  of  the  state  of  this  knowledge. 

The  probable  pilot  behavior  for  any  given  vehicle  flight  condition  and 
task  prior  to  the  PIO  can  be  estimated  using  the  model  of  Fig.  2  directly, 
although  some  simplifications  are  usually  made  (see  Refs.  9-1 6,  and 
especially  Refs.  13  and  1 6) •  The  essence  of  the  simplified  adaptation 
rulea  is  that  the  pilot  adopts  a  form  of  equalization  such  that  the  pilot- 
plus-vehicle  open-loop  describing  function  is  much  greater  than  1 .0  at  low 
frequencies  while  approximating  a  -20  db/decade  slope  in  the  region  of 
crossover  frequency,  cq,.  These  adjustments  satisfy  the  conflicting  demands 
of  good  closed-loop  control  at  the  input  frequencies  lower  than  cq>  (below 
afc,  |yol|  »  1$  thus  | YqlI  =  1.0),  and  adequate  stability  (at  cq.  an  adequate 
phase  margin  exists).  Recent  experiments  by  Systems  Technology,  Inc., 
and  The  Franklin  Institute  (Ref.  25)  show  that  under  laboratory  conditions, 
the  crossover  frequencies  for  a  given  controlled  element  are  constant  so 
long  as  tq,  >  (where  is  the  effective  bandwidth  of  the  system  forcing 
function).  This  constancy  of  iq»  does  not  extend  across  different  controlled 
elements,  but  ranges  from  3  to  8  rad/sec,  the  lower  for  control  of  pure 
inertias  and  the  higher  for  control  of  pure  gain  elements.  Comparison  of 


measured  in  flight  versus  simulator  shows  that  about  half  the  laboratory 
value  of  at  is  realized  (Ref.  58).  In  general,  the  simplified  rule  of  thumb 
for  pre-PIO  adaptation  is:  Make  |YqjJ  *  | Yp  •  Yci  =  | a^/ s |  in  the  vicinity 
of  at  «  1 .0  to  2.0  rad/sec. 

For  the  second,  or  actual  PIO  phase,  the  dominant  loops  must  be  treated 
with  sinusoidal-input  describing  function  models.  The  state  of  knowledge 
for  this  periodic-input  case  is  almost  the  opposite  of  that  for  random 
inputs.  Here  a  great  deal  is  known  about  sinusoidal- input  describing  func¬ 
tions  for  inanimate  nonlinear  elements,  but  some  critical  information  is 
lacking  on  sinusoidal-input  models  for  the  pilot.  The  present  hypothetical 
pilot  describing  function  model  for  PIO  derives  from  sinusoidal- input 
frequency  response  measurements  in  manual  control  systems  (Refs.  26-31 ) . 

In  these  studies,  conducted  with  a  pure  gain  controlled  element,  a  key 
observation  is  that  the  terminal  phase  of  pilot  adaptation  is  "synchronous," 
or  "precognitive"  behavior.  In  other  words,  once  the  pilot  recognizes  that 
his  input  is  a  sinusoid,  he  essentially  duplicates  the  sinusoid  with  no 
phase  lag — basically  Yp  =  Kp.  This  nondelayed  behavior  breaks  down  above 
3  cps  or  so. 

PIC  investigations  using  a  Yp  »  Kp,  with  visual  inputs  (e.g,,  pitch 
angle)  presumed,  have  had  a  distinct  measure  of  success  (e.g..  Ref.  1,  7> 
and  20)  in  determining  probable  causes  and  assessing  cures  for  existing 
PIO’s,  Despite  these  successes  the  confidence  level  in  the  Yp  =  Kp  form 
for  the  human  pilot  sinusoidal-input  describing  function  is  not  too  high. 

For  instance,  in  most  PIO's  the  pilot  has  adapted  to  the  pre-PIO  situation 
with  lead  or  lag  or  lag-lead  equalization  generated  as  the  occasion  demands. 
The  question  then  arises  as  to  whether  the  pilot  retains  the  equalization 
as  the  PIO  begins,  or  whether  he  abandons  the  equalization  and  attempts 
synchronous  behavior  with  essentially  zero  phase  lag.  The  effective 
reaction  time  delay  no  longer  exists  in  his  response.  Also,  for  a  pure 
gain  controlled  element  the  usual  low  frequency  lag  equalization  present 
for  random  inputs  is  dropped  for  the  sinusoidal- input  case.  But  it  is 
not  clear  whether  a  net  phase  angle,  made  up  of  the  sum  of  the  leads  and 
lags  in  the  previous  equalization,  is  dropped  or  retained  in  the  PIO  for 
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controlled  elements  other  than  Yc  *  Kq.  This  question  is  important  in 
determining  the  gain  reduction  required  to  re-establish  stable  control 
and  in  estimating  whether  or  not  a  PIO  is  likely  to  occur  for  a  specific  * 
flight  situation.  New  experimental  research  is  needed  on  this  matter. 

Between  the  pre-PIO  and  developed-PIO  phases  is  a  transient  period 
which  contains  the  excitation,  vehicle  change,  pilot-organized  system 
change,  etc.,  which  serves  to  initiate  the  PIO.  This  transient  takes  a 
finite  time,  during  which  all  stationary  descriptions  of  pilot  behavior 
are  invalid.  The  only  method  which  can  be  used  at  this  state  of  the  art 
is  to  analyze  the  closed-loop  dynamics  both  before  and  during  a  PIO  to 
see  if  conditions  conducive  to  sustained  oscillations  exist. 

D.  PILOT  INPUTS  AND  OUTPUTS  FOR  PIO 

The  easiest  way  to  represent  the  pilot  inputs  and  outputs  is  on  a 
block  diagram.  A  typical  block  diagram  of  a  pilot-vehicle  system  is 
shown  in  Fig.  3  for  longitudinal  PIO,  and  similar  block  diagrams  can  be 
drawn  for  other  axes.  Representative  mathematical  models  for  various 
system  elements  for  a  fighter-type  aircraft  are  shown  at  the  bottom  of 
Fig.  3. 

The  inputs  to  the  pilot  prior  to  a  PIO  consist  of  those  visual  and 
motion  cues  which  are  used  to  perform  the  basic  flying  task  at  hand. 

Most  of  the  complex  piloting  tasks  involve  aircraft  attitude  relative  to 
a  real  or  artificial  horizon  as  the  inner  loop.  In  other  cases,  such  as 
climb  programming  or  LABS  maneuvers,  the  normal  acceleration  indicator  may 
be  in  dominant  use.  Another  task  in  which  PIO  is  frequently  encountered 

is  formation  flying,  where  the  cues  are  a  complex  combination  of  linear 
displacements,  attitudes,  and  angular  rates. 

Once  a  PIO  limit  cycle  has  been  established,  the  input  situation  is 
not  so  well  understood  experimentally.  However,  one  of  the  dominant 
inputs  is  certainly  the  physical  acceleration  felt  by  the  pilot,  and 
others  are  the  gross  motions  of  the  outside-world  horizon.  Experience 
has  shown  that  most  of  the  PIO  situations  and  pilot  behavior  can  be 
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understood  by  considering  the  attitude  and/or  normal  accelerations  as 
the  dominant  sensed  quantities. 

The  pilot  outputs  of  interest  axe  the  forces  and  displacements  applied 
to  his  "manipulator,"  i.e.,  control  stick  or  pedals.  As  Fig.  3  shows,  a 
pilot  has  both  force  and  displacement  control  loops  within  his  neuro¬ 
muscular  actuator  system.  A  great  deal  of  (largely  intuitive)  controversy 
has  raged  over  whether  the  pilot  output  should  be  considered  to  be  essen¬ 
tially  a  force  or  essentially  a  displacement.  In  the  more  sophisticated 
context  of  Fig.  3,  this  resolves  to  the  question  of  the  relative  degrees 
of  "tightness"  of  the  force  and  displacement  loops.  An  answer,  with  its 
qualifying  details,  has  not  yet  been  obtained,  although  some  experimental 
programs  are  under  way  to  help  establish  the , fundamental  neuromuscular 
processes  involved.  The  pertinent  main  points  which  have  evolved  from 
extensive  experience  to  date  are: 

1 .  Both  force  (kinesthetic)  and  displacement  (proprio¬ 
ceptive)  feedbacks  are  used  for  control  by  the 
neuromuscular  system. 

2.  With  spring-restrained  manipulators  the  displacement 
outputs  seem  to  be  dominant  (high  displacement  loop 
gain)  when  large  or  preprogrammed  manipulator  motions 
are  required,  whereas  the  force  outputs  are  dominant 
(low  displacement  loop  gain)  when  fine  motions  and 

,  precision  control  are  needed. 

3.  Even  with  tight  closure  of  the  displacement  loop, 
such  as  occurs  when  a  free-moving  pure  displacement 
manipulator,  the  high  frequency  dynamic  characteris¬ 
tics  of  the  neuromuscular  system  show  more  lag  than 
those  evident  with  a  "stiff-stick"  pressure  control. 

•The  second  point  above  qualitatively  agrees  with  experience  on 
longitudinal  and  lateral  controls.  Thus  the  force-feel  characteristics 
of  the  aircraft  control  system  seem  to  be  more  critical  than  the  displace¬ 
ment  characteristics  for  small-motion  longitudinal  control.  This  implies 
that  the  pilot's  neuromuscular  system  displacement  loop  is  operating  at 
low  gains  since,  otherwise,  a  tight  proprioceptive  loop  would  suppress 
many  of  the  control  system's  objectionable  force  nonlinearities.  In 


aileron  control,  on  the  other  hand,  control  motions  can  be  more  gross 
and  less  emphasis  on  the  minimization  of  force-type  nonlinearities  is 
needed. 

Summarizing  these  results,  it  appears  that  either  force  or  position 
outputs  may  be  dominant  in  PIO  studies,  depending  on  the  situation  being 
considered  and  the  nonlinearities  present..  One  should,  therefore,  check 
the  pilot-vehicle  system  characteristics  for  both  force  and  displacement 
inputs  from  the  pilot  to  see  if  conditions  for  a  sustained  limit  cycle 
can  be  found.  Further  experiments  related  to  these  problems  are  planned. 
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Certain  necessary  steps  in  a  sequence  are  common  to  most  PIO's.  Before 
discussing  these  from  the  standpoint  of  PIO  causes,  consider  first  some  of 
their  physical  attributes.  Figure  1  clearly  illustrated  the  following 
sequence : 

1 .  The  pilot  is  initially  adapted  to  some  control 
situation  between  the  extremes  of  precise  tracking 
maneuvers  and  hands-off  flight. 

2.  Something  causes  a  sudden  change  in  the  situation 
dynamics.  This  could  be  a  change  in  the  pilot's 
organization  of  the  system,  the  initiation  of  a 
large  steady  maneuver,  a  damper  failure,  etc.  (in 
Fig.  1 . the  change  followed  the  shutdown  transient 
of  a  faulty  pitch  damper,  which  left  the  elevator 
in  a  mistrimmed  position.) 

3.  An  oscillation  builds  up  and  is  sustained  for  a  few 
cycles ;  a  legitimate  limit  cycle  exists.  (Note 
that  two  sets  of  limit  cycles  exist  in  Fig.  1  for 
times  greater  and  less  than  15  sec.) 

4.  Finally,  the  pilot  either  lets  go,  freezes  the 
stick,  or  puts  in  a  well-timed  maneuver  to  kill  the 
oscillation. 

The  conditions  that  lead  to  step  3  that  are  the  root  "cause"  of  a  PIO. 
It  is  always  possible,  by  making  the  forward  loop  gain  high  enough,  to 
drive  any  physical  closed-loop  system  unstable.  Pilot  trainees  or  test 
pilots  feeling  out  a  new  aircraft  often  tighten  up  on  their  control 
response  enough  to  provide  one  or  two  oscillations  indicating  incipient 
instability.  If  this  tendency  is  easy  to  avoid,  and  if  a  modest  reduction 
in  pilot  gain  (e.g.,  25  to  50  percent)  can  remove  the  instability,  then 
such  situations  do  not  usually  end  up  as  serious  PIO  cases.  It  is  those 
unusual  situations  where  several  factors  combine  to  make  the  pilot— vehicle 


instability  region  either  impossible  to  avoid  or  difficult  to  get  out  of 
that  are  the  real  concern  of  this  report. 

There  are  several  instances  where  a  sustained  oscillation  of  the 
vehicle-plus-control-systems  alone  may  be  triggered  by  a  large  pilot- 
induced  maneuver  (e.g.,  bobweight  plus  nonlinear  gearing).  Although  such 
oscillations  do  not  involve  visual  or  acceleration  feedbacks  through  the 
pilot's  higher  centers,  they  often  depend  on  the  presence  of  the  pilot's 
am  on  the  stick.  Thus  the  pilot's  neuromuscular  system  is  involved,  and 
the  oscillations  are  still  "caused"  by  the  pilot.  Such  oscillations  are 
often  unavoidable  because  it  is  seldom  practical  to  completely  let  go  of 
the  stick,  especially  at  the  high- subsonic,  low-altitude  conditions  where 
the  problem  is  aggravated. 

B.  PXLOT-REJERHICED  CAUSES  OF  PIO 


There  are  several  ways  of  looking  at  the  causes  of  a  PIO.  One  is  to 
catalog  all  the  PIO  situations  ever  recorded,  Including  all  the  necessary 
,  subsystem  details,  etc.,  and  then  to  say  that  each  combination  of  vehicle 

and  subsystems  when  combined  with  the  pilot  was  the  cause  of  a  PIO. 

Another  way  is  to  note  that  certain  system  phenomena,  such  as  stick-force— 
to-control-def lection  hysteresis,  often  lead  to  PIO  when  other  conditions 
are  right  and  can  thus  cause  PIO.  A  third  way,  and  one  which  seems  to 
transcend  the  difficulties  of  the  previous  two,  is  to  say  that  certain 
inherent  human  physical  limitations  are  the  basic  cause  for  any  PIO.  This 

I 

is  not  to  degrade  the  human  pilot's  role  but,  instead,  to  emphasize  it, 
because  it  is  unlikely  that  any  black-box  could  be  devised  which  is  as 
clever  and  effective  in  coping  with  unmanageable  controlled  elements  as  a 
skilled  pilot.  Were  it  not  for  the  pilot's  versatile  gain  adaptability, 
many  flight  conditions  would  be  unstable.  But  there  is  a  limit  to  the 
rapidity  with  which  the  human  can  adapt,  and  thi  s  can  sometimes  lead  to 
a  PIO. 

When  referred  to  the  pilot,  then,  the  basic  causes  of  PIO  seem  to  fall 
into  the  following  citegories: 

( 
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Incomplete  pilot  equalization 

a.  Incomplete  training 

b.  ■  Inappropriate 'transfer  of  adaptation 
(i.e.,  carryover  of  improper  techniques 
from  another  aircraft) 

♦ 

Excessive  demands  on  pilot  adaptation 

a.  Required  gain,  lead,  or  lag  lie  outside 
the  range  of  normal  capabilities 

b.  Rate  of  adaptation  is  too  slow  to  preclude 
oscillation 

c.  Inadequate  capability  to  cope  with  system 
nonlinearities 

3«-  Limb-manipulator  coupling 

a.  Impedance  of  neuromuscular  system  (including 
limb)  on  control  stick  or  pedals  changes 
feel  system  dynamics 

b.  Motion-induced  limb  force  feedback  (e.g. , 
arm  becomes  a  bobweight) 

Incomplete  Pilot  Equalization  is  a  common  problem  and  results  from 
insufficient  time  or  trials  with  the  given  situation  for  the  pilot  to 
achieve  a  good,  stable  closure.  The  cure,  of  course,  is  more  training, 
and  more  than  one  PIO  problem  has  disappeared  as  the  piloting  technique 
is  refined.  In  sor.e  instances  the  required  piloting  technique  may  be 
quite  different  from  that  with  which  the  pilot  is  familiar,  and  a  diffi¬ 
cult  transition  period  results.  If  the  analyses  show  that  no  equalization 
is  needed  by  the  pilot  prior  to  a  PIO,  but  that  the  required  pilot  gain  is 
very  low  (i.e.,  control  is  very  sensitive),  then  incomplete  pilot  gain 
adjustment  may  well  cause  PIO  tendencies.  Due  to  the  temporary  nature  of 
these  tendencies,  this  cause  of  PIO  is  not  considered  as  serious  as  some 
of  the  others. 

Excessive  Demands  on  Pilot  Adaptation  are  the  most  common  cause  of 
malignant  PIO's.  The  ranges  of  lead,  lag,  and  gain  available  to  the 
pilot  are  fairly  broad,  out  tiiere  are  some  limitations  which  can  result 
in  a  sustained  oscillation.  Experience  on  fixed-base  flight  simulators 
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indicates  that  a  pilot  can  change  his  gain  over  nearly  two  decades  if 
given  enough  time  (Ref.  13)*  However,  for  a  sudden  change  in  the  required 
pilot  gain,  a  2:1  variation  (±6  db)  is  easily  accomplished,  but  a  4:1 
change  (±2  db)  is  apparently  quite  difficult,  at  least  for  aircraft 
control  tasks.  The  achievable  lead  time  constants  measured  in  the  STI-FIL 
experiments  range  from  Tj,  =  0  to  greater  than  3  sec,  while  achievable 
lags  range  from  Tj  =  0  to  10  sec.  The  pilot  can  easily  adopt  small  leads 
roughly  equal  to  his  reaction  time  delay,  but  the  generation  of  small  lags 
in  the  vicinity  of  0.1  to  1  sec  requires  relatively  longer  training.  As 
noted  previously,  when  a  transient  oscillation  starts  to  dominate  the 
input,  the  pilot  seems  to  reverb  to  a  pure  gain  behavior,  thus  dropping 
his  previous  equalization.  In  certain  situations,  where  the  pilot  has 
adopted  a  lag  equalization  to  compensate  for  a  high-frequency,  low-damped 
pure  second-order  controlled  element,  dropping  the  lag  equalization  can 
seriously  destabilize  the  system,  leading  to  a  series  of  quasi-steady 
PIO's  as  the  pilot  is  forced  to  backpedal  severalfold  on  his  gain  to 
restabilize  the  system.  This  cause  of  PIO  may  be  difficult  to  prove  at 
the  present  state  of  the  art. 

By  far  the  most  excessive  demands  on  pilot  adaptation  are  caused  by 
various  nonlinearities  in  the  control  system.  Certain  forms  of  complex 
series  nonlinearity,  such  as  elevator-to-stick-force  hysteresis  or  rate 
limiting  Inside  a  closed-loop  system,  defy  the  pilot’s  ability  to  compen¬ 
sate  or  invert  it,  and  some  degree  of  PIO  becomes  inevitable  during  tight 
closed- loop  control.  These  problems  can  be  readily  treated  and  understood 
by  analytical  methods  to  be  described  later. 

The  last  category  of  PIO  causes  is  related  to  Limb-Manipulator  Coupling 
This  is  a  general  term  for  certain  dynamic  interactions  such  as  the  effect 
of  the  pilot's  arm  mass  acting  like  a  bobveight  and  thus  feeding  back  local 
accelerations  into  the  control  system.  This  type  of  coupling  can  actually 
destabilize  the  short-period  mode  if  the  feel  system  natural  frequency  is 
nearby,  and  if  the  pilot  does  not  fight  the  acceleration  loads  or.  his  arm. 

A  more  subtle,  but  possible,  PIO  cause  is  the  mass-like  impedance  o’f  the 
pilot’s  arm  (or  legs)  on  the  controls  when  he  is  loosely  hanging  on  to  the 
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i  }  stick.  Depending  on  the  apparent  arm  mass  "seen"  by  the  feel  system,  its 

natural  frequency  can  be  lowered  enough  that  coupling  forces  (e.g.,  bob- 
weight)  can  drive  the  short-period  feel  system  combination  unstable 
(Refs.  1  and  7).  Both  types  of  limb-manipulator  effects  need  further 
investigation  using  moving-base  simulators. 

In  the  next  section  it  is  shown  thal^  in  classifying  the  types  of 
physical  phenomena  and  analysis  techniques  associated  with  PIO,  there 
appears  to  be  a  natural  division  into  linear  and  nonlinear  categories. 
Correlating  the  examples  given  later  in  Table  I  with  the  basic  causes 
just  described,  reveals  that  the  "linear"  types  of  PIO  are  usually  caused 
by  the  first  category  (incomplete  Pilot  Equalization),  while  the  nonlinear 
types  of  PIO  ell  fall  into  the  second  category  (Excessive  Demands  on  Pilot 
Adaptation) .  Although  the  correspondence  should  not  be  drawn  too  closely, 
it  is  not  fortuitous  and  was  one  reason  for  the  scheme  chosen  to  categorize 
the  types  and  causes  of  PIO. 

('  ) 
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SECTION  IV 
AX4LYBIS  07  PIO 


A.  GENERAL  APPROACH 

The  object  of  PIO  analysis  is  to  determine  the  fundamental  controlled 
element  problems  which  could  cause  an  oscillation  when  the  pilot,  repre¬ 
sented  by  a  describing  function,  attempts  to  close  the  control  loop.  The 
basic  approach  to  the  analysis  of  PIO  is  to  select  a  potentially  critical 
aircraft  and  flight  condition,*  then  to  analyze  this  condition,  using 
servoanalysis  techniques  and  various  pilot  models  representing,  reasonable 
piloting  behaviors,  to  determine  whether  either  zero  damping  (in  a  linear 
system)  or  a  limit  cycle  (for  a  nonlinear  system)  can  exist.  When  a 
diagnostic  analysis  is  required,  the  criticax  flight  conditions  where  PIO 
are  suspected  or  experienced  are  usually  given  at  the  outset,  although 
other  flight  conditions,  with  no  x'ecord  of  PIO  (or  in  regions  where  no  PIO 
is  expected)  should  also  be  examined  to  check  the  validity  of  the  results. 
For  a  prognostic  (at  the  design  stage)  analysis,  election  of  potential 
PlO-prone  flight  and  control  system  conditions  is  primarily  a  matter  of 
identifying  situations  which  contain  one  or  more  characteristics  inimical 
to  pilot-vehicle  system  stability.  To  some  extent  this  can  be  made 
deductive,  e.g.,  by  cataloging  conditions  representing  dynamic  response 
extremes  in  the  control  system  and  vehicle;  but  experience  and  knowledge 
of  past  PIO  history,  with  its  indications  of  possible  parallels,  provide 
the  best  guide.  This  will  be  discussed  further  in  the  next  subsection. 

The  analyst  is  cautioned  to  examine  all  possibilities  for  PIO  before 
settling  on  the  dominant  cause,  since  experience  has  shown  that  more  than 
one,  and  usually  a  combination  of,  PlO-enhancing  factors  are  involved  in 
say  particular  problem. 

In  subsections  to  follow  PIO's  are  classified,  and  the  corresponding 
critical  systems  and  conditions  are  discussed.  Next,  some  typical  steps 
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[  and  data  required  in  any  PIO  analysis  are  reviewed.  This  material  will 

help  to  define  the  scope  of  test  and  analytical  work  required.  Finally, 
four  examples  of  distinctly  different  PIO  problems  will  be  presented  to 
illustrate  the  range  of  problems,  the  techniques,  and  the  forms  of 
mathematical  models  which  are  involved  in  PIO  analyses. 

B.  CLA00Z7ICASZON  07  PIO 

The  profusion  of  individual  cases  of  PIO  which  have  been  and  will  be 
encountered  requires  some  unifying  classification  scheme.  Of  the  several 
alternatives  possible,  the  following  scheme  seems  to  be  the  best  suited 
from  the  point  of  view  of  the  physical  phenomena  involved  and  the  type  of 
analysis  required: 

PIO's  are  classed  by  the  control  axis  involved,  i.e,,  as 
longitudinal  (pitch) ,  lateral  (roll),  directional  (yaw), 
and  coupled  lateral-directional. 

The  type  is  defined  by  the  nature  of  the  physical  phenomena 
(and,  incidentally,  the  kind  of  analysis)  involved: 

>r 

^  I.  Oscillations  due  to  linear  pilot-vehicle 

coupling 

II.  Limit  cycles  due  to  oue  or  more  nonlinear 
elements  in  series  in  the  primary  control 
loop 

III.  Limit  cycles  due  to  one  or  more  nonlinear 
elements  in  vehicle  motion  feedback  paths 
subsidiary  to  the  primary  conti’ol  loop 

The  species  of  PIO  within  each  class  and  type  describe 
the  specific  character  of  each  individual  case  of  PIO. 

From  an  analysis  standpoint,  the  class  distinction  directly  indicates 
the  appropriate  equations  of  vehicle  motion.  The  type  iefinition  orders 
the  problem  in  terms  of  analytical  rather  than  physical  complexity.  Thus 
Type  I  systems  require  only  linear  analysis  methods,  although  the  systems 
may  be  multiloop  in  nature,  whereas  Type  II  and  Type  III  systems  always 
require  nonlinear  analysis  techniques.  Type  II  systems  are  essentially 
single-loop  insofar  as  the  nonlinear  analysis  is  concerned,  but  Type  III 
systems  are  always  multiloop.  The  species  description  is  intended  to 
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emphasize  physical  rather  than  analytical  characteristics,  so* it  may 
provide  only  anecdotal  inf o.rmation* 

As  an  example  of  the  classification  scheme,  the  "J.  C.  Maneuver" 
would  be  categorized  as  an  oscillation  in  the  pitch  axis  (class)  involving 
nonlinearities  in  the  feel  system  (type),  and  characterized  by  a-rapid- 
ospillation-at-large**ioad-factors-during-flight-maneuvers-at-high-dynamic- 
pressure  (species).  Several  known  or  suspected  PIO  cases  are  classified 
using  the  above  basesrin  Table  I.  As  a  partial  substitute  for  direct 
experience,  Table  I  can  be  of  significant  help  with  either  diagnostic  or 
prognostic  PIO  analyses.  In  the  former  case,  when  the  nature  of  the  PIO 

can  be  identified  with  one  of  those  shown,  the  available  references  and 

/ 

type  of  analysis  required  are  immediately  apparent.  In  the  latter  case, 
when  predictions  of  potential.  PIO  problems  for  a  new  design  are  required, 
Table  I  shows  the  most  likely  possibilities,  as  well  as  those  amenable  to 
simple  analysis  at  the  preliminary  design  stage.  Also  shown  in  Table  I 
are  the  critical  flight  conditions  and  critical  subsystems  for  each  of  the 

species  of  PIO  listed.  Each  listing  follows  the  order  Species  (Air- 

Ref  ”* 

craft)  * }  Critical  Subsystem;  Critical  Flight  Condition;  Remarks,  and 

makes  use  of  the  shorthand  legend  identified  at  the  bottom  of  the  table. 

0.  DAS A  AND  8TE78  NWXXED 

Experience  has  shown  that  the  detailed  nature  of  the  subsystem  dynamics 
must  be  at  hand  for  a  meaningful  analysis  to  be  undertaken.  For  the  simplest 
type  of  PIO  (those  involving  linear  pilot-vehicle  coupling)  only  one  more- 
or-less  unique  transfer  function  for  each  subsystem  need  be  known.  However, 
many  PIO  involve  subsystem  nonlinearities  of  varying  degrees  of  complexity. 
For  such  instances  sinusoidal- input  describing  function  techniques  have 
proven  to  be  most  useful,  and  so  it  becomes  essential  to  obtain  the  various 
subsystem  describing  functions  for  various  levels  of  input  to  the  nonlinear 
elements.  Table  II  has  been  prepared  as  a  guide  to  the  scope  of  information 
required  for  each  critical  subsystem  as  indicated  in  Table  I.  The  items  in 
Table  II  are  by  no  means  exhaustive,  and  serve  merely  to  indicate  the  scope 
of  information  required  and  the  form  in  which  it  will  prove  most  useful. 
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TABLE  I 

CLASSIFICATION  OF  SOKE  KNOWN  710  CASES 

Examples  .hdwn.aat  SPECIES  (Aircraft)2'^*;  Critical  Subsystem;  Critical  Flight  Coalition;  Peoria 


CLASS 

TYPE 

I.  .  LINEAR 

FITCH 

IMPROPER  SIMUIATION^j  D,  Vj  a:  Abnormally  high  value  of  1/T(?2  end  lov  Cc%p  led  to  aero  gBp  when 

regulating  large  disturbances. 

GCA-INDUCED  PHDGOID  (C-97)*^i  D,  V;  b:  lag  from  radar-do  tec  ted  error  to  vol<*»  <*Ai»»nd  led  to 

unetable  closed-loop  phugold  node, 

ARM  ON  STICK  (AhD-1 )^.  (T-3&A)'  i  Ft  at  Arm  mass  Increases  feel  system  inertia;  leads  via  B  feed- 

hack  to  unstable  coupling  with  short-pi.rlod  dynamics  if  pilot  merely  hangs  loosely  onto  stick 
after  a  large  input. 

IATERAIrDIRECTIONAL 

V“H  EFfECT  (X-15)20,  (T-3JVSA)2’,  (F-101B;J\  (F-tO&),  (KC-135A)55,  (B-56)i  V;  ct  Zeros  of 

roll/aileron  transfer  function  are  higher  than  Dutch  roll  frequency,  |iqj/<qj|  >  i.O,  leading  to 
..closed-loop  instability  at  lov  conditions. 

BORESIGHT  OSCILLATIONS  (F^tjA )i  D,  V;  ct  Spiral  roll  mode  driven  unstable  if  roll  infonmtion  is 
degraded  during  gunnery. 

YAW 

FUEL  SL06H  SNAKING  (KC-IKA)55.  (T-BfA)52;  V;  ct  Fuel  aloah  mode  couDle.  vlth  Dutch  roll  mode 

when  rudder  used  to  atop  yaw  oscillation. 

ROLL 

NONE  KNOWN 

Legend; 


Critical  Subsystems t 


Critical  Flight  Conditional  • 


Superacrlpta  refer  to 
reference!  of  occurrence. 


D  ■  Display  a 

F  •  Feel  system  (except  B)  h 

B  x  Bobvelght 

S  ■  Power  aervo  actuator  c 

V  ■  Vehicle  (airframe) 

A  ■  Augmenter  (damper; 


»  Low  altitude,  near-!onlc  Mach 
»  Landing  approach  or  takeoff 
■  Crulee 
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CLASS 


TYPE 


II.  SERIES  NONLINEAR  ELEMENTS 


PORPOISING  (SB2C-1)^;  Pi  ci  Hyateresii  in  stick  vsrsus  slsvster  deflection  resulted  in  lov 
frequency  speed  end  climb  sscUIstlons. 

J.  C.  MANEUVER  (F-B6D),  (r-IOOC)i  T,  Si  s)  Velv*  friction  plus  cenplient  cabling  resulted  in 
large  oscillations  at  abort  period. 

PITCH-UP  fxy-IQb,  (r-IOlB)3**.  (P-102A)^;  Vi  cl  Unstable  kink  in  M(a)  curve  led  to  moderate-period 
oscillations  of  varying  amplitudes  (depending  on  extent  and  nature  of  the  kink)  during 
aaneuvers  near  the  crtlcal  angle  of  attack. 


bt  Closed-loop  around  elevator  rate-limiting  cauaed  aoderate  oscilla¬ 


tions  at  short  period. 


UTOAU-DIREmOk'U, 


TRANSONIC  SNAKING  (A5D)l  V,  Pi  a,  ci  Separation  over  rudder  causes  control  reversal  for  snail 
deflections,  leading  to  Unit  cycle  if  rudder  used  to  daop  yav  oscillations. 


PILOT- INDUCED  CHATTER  (F-10bB 
attempted  to  control  It. 


Ai  cl  Snail  Unit  cycle  due  to  damper  aggravated  vhenever  pilot 


III.  SUBSIDIARY  EEEDBACX  NONLINEAR  ELEMENTS 


WEI0HT  BREAKOUT  (ALP-1  (T-38A)1,lfi  T,  Bi  ai  At  high-g  maneuver#  the  bolvelght  overcomeo 
eyatea  friction  and  reduces  apparent  damping  of  the  aircraft  in  response  to  force  Inputs, 
resulting  in  large  oscillations  at  abort  period. 


LOSS  OF  PITCH  DAMPER 


LOSS  OP  YAW  DAMPER 


IATERAIHJIRECTIONAL 
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TABLE  II 

SUBSYSTEM  DATA  REQUIRED  FOR  PIO  ANALYSES 


SUBSYSTEM  DATA  REQUIRED 

D  =  Displays  Type  of  presentation:  Compensatory  versus  pursuit. 
Instrument  dynamics:  6-q/Q,  n^/n^,  etc. 

Kinematics  of  visual  field*:  e.g.,  horizon  visibility, 
orientation  in  formation  flight,  visual  landing  aids. 


F  =  Feel 
System 
(excl.  B) 


Control  surface  commands  due  to  stick  force  and 
displacement : 

Frequency  response,  6C/FS  =  f|Fs|,  5C/6S  =  f|6s|; 
three  levels  of  Fs  or  os,  for  at  least  one  decade 
centered  at  o^p. 

Static  data  on  friction,  hysteresis  loop  in  5C  versus 
Fs,  etc. 


For  nonlinear  stick  gearing,  take  describing  function 
data  near  &ctriln* 


B  =  Bobweight 


Bobweight  force  feedback  dynamics:  Fs  /ai,  Fs  /&*,  etc. 

Bobweight  breakout  and  hysteresis  as  f(n£),  etc. 

Equivalent  bobweight  location,  considering  all  mass 
unbalance  which  is  sensitive  to  nz  and 

Possibly  include  inertial  bending  and  support  deflec¬ 
tion  under  nz  leads  as  B  effects,  if  not  already  in 
vehicle  dynamics. 


S  =  Servo 
Actuator 


Control  output/input  frequency  response  and  describing 
function: 


8/8c  at  three  levels  of  5C;  also  6e  ,  , 
threshold,  etc. 


-max 


> 


A  =  Augmenter 
(Damper) 

i 


Control  commands  due  to  vehicle  motions: 

5c^/body  rates;  describing  functions  if  nonlinear; 
threshold,  limits,  etc. 


V  -  Vehicle 
(Airframe) 


Specific  configuration  and  flight  conditions: 

Compatible  weights,  inertias,  a,  q,  5-trim>  ^trim' 
Mtrim;  etc. 


Calculate  transfer  functions  or  measure  frequency 
response  (for  constant  | 5 | ) : 

a^/8  (aj  =  accel.  at  c.g.,  pilot,  effective  B  loca¬ 
tion,  etc,);  attitude/6  (attitude  =  0  or  <p  or  \|f,  etc., 
or  corresponding  rates). 

Preferably  include  inertial  and  aeroelastic  bending  here; 
also  Zge  or  Y5  effects. 


2k 


The  principal  steps  involved  in  any  analysis  of  pilot-induced 
oscillations  are  as  follows: 


1.  Select  the  dominant  control  node  (e.g.,  pitch 
attitude  control) ,  critical  flight  conditions, 
and  dominant  subsystems  likely  to  be  directly 
involved  in  any  oscillation  (see  Tables  I  and  II). 

2.  Obtain,  either  by  analysis  of  measurement,  the 
required  transfer  functions  or  describing'  func¬ 
tions  of  the  critical  systems  as  specified  in 
Table  II. 

3*  Assume  an  oscillating  output  at  selected  levels 
(either  observed  or  estimated),  then  perform  a 
closed-loop  describing  function  analysis  to 
verify  if  the  conditions  required  to  sustain  the 
osr. illation  can  be  reasonably  expected.  Two  or 
more  levels  of  oscillation  may  be  necessary  to 
cover  the  extremes  of  the  subsystem  nonlinear 
describing  functions.  Use  a  pure  gain  model  for 
the  pilot. 

4.  Check  to  see  whether  the  pilot  gain  required  in 
the  limit  cycle  is  consistent  with  the  probable 
magnitude  of  the  pilot’s  Gaussian-input  describing 
function  adapted  prior  to  the  onset  of  PIO  (magni¬ 
tude  of  Yp  at  the  PIO  frequency).  The  methods  of 
Ref.  13  or  16  should  be  used  to  estimate  the 
pre-PIO  pilot  adaptation  for  the  specific  control 
task  involved. 

5.  If  conditions  for  an  oscillation  are  compatible 
with  Steps  3  and  4,  then  a  PIO  is  likely,  provided 
the  assumed  mode  of  pilot  control  (e.g.,  attitude 
error  to  force  output)  is  valid. 

6.  The  analysis  should  be  duplicated  for  any  other 
likely  pilot  control  modes  or  nonlinearity  ampli¬ 
tudes,  and  the  predicted  PIO  tendencies  weighted 
according  to  the  results  therefrom. 


D.  EXAMPLE  CASKS  , 

The  examples  to  follow  will  illustrate  potential  PIO  situations  of 
increasing  analytical  complexity,  corresponding  to  the  type  of  classifi¬ 
cation  given  above.  The  first  two  cases  invoV^  no  nonlinearities  what¬ 
ever,  so  both  illustrate  Type  I  situations.  The  third  example  considers 
a  piloted-force  nonlinearity  which  is  analytically  tractable;  the  non¬ 
linear  system  element  is  also  in  series  with  the  primary  control  and  hence 
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the  control  situation  is  Type  Hi  Finally,  the  fourth  example  treats  a 
TJrpe  III  system  which  is  so  complex  that  it  can  be  feasibly  analyzed  only 
at  the  two  extremes  of  nonlinear  behavior. 

1 .  Erjuqplt  1 :  Linear  Control  of  Short-Period  Longitudinal  Motions 

Consider  the  pitch  control  of  an  aircraft  at  high  subsonic  speeds 
and  lew  altitude.  As  noted  earlier,  in  a  sustained  PIO  the  pilot's 
sinusoidal  describing  function  can  be  considered  a  simple  gain  (Yp  *  Kp) . 
Under  such  circumstances  end  assuming  that  hi6  primary  response  is  to 
1  visual  pitch-attitude  cues,*  the  pilot-vehicle  open-loop  describing  func¬ 
tion  at  any  limit  cycling  frequency,  o,  is  given  by 


where  s  «  jn  and  phugoid  motions  are  neglected.  Accordingly,  the  closed- 
loop  characteristic  equation  is  given  by 

A'  -  s|s2  +  2£SpU)5pS  +  o|p|  +  RpM^s  + 

-  +  sa(2{Sp<%p)  +  “(<4  +  Kp*>&)  +  (2) 

The  usual  factored  form  of  A'  is  expressed  as  follows,  where  the  prime 
denotes  closed-loop  parameters: 

*It  should  be  emphasized  here  that  in  spite  of  accompanying  vertical 
accelerations,  attitude  cues  will  be  those  the  pilot  consciously  uses  in 
attempting  to  get  out  of  a  PIO  situation,  so  the  above  considerations  are 
valid.  Nevertheless,  acceleration  inputs  introduced  by  the  pilot's  arm 
bobweight  effect  acting  through  control  system  friction  could  be  a  non¬ 
linear  destabilizing  influence.  The  fourth  example  considers  such  effects 
in  detail. 
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A* 


{,+k)<f*H  »pV  +  “»p) 

(i  +  2?'  \  +  sL'2  +  £&2k)  + 

\TC  ^sp^sp  /  \“6P  Tc  /  Tc 


=  SJ  +  s 


(3) 


Since  a  sustained  oscillation  implies  zero  damping,  the  conditions  for 
PIO  (and  the  only  condition  for  which  Eq  3  has  meaning,  because  of  the 
assumed  oscillation)  correspond  to  those  for  which  £spc%p  is  zero. 
Equating  the  s2  coefficients  of  Eq  2  and  3> 


2^spa)sp  +  T„  “  2^sp% 


and  the  condition  for  an  oscillation  becomes 


25spcusp  "  2^Bpmsp  tc  *  0  M 

• 

Thus,  if  l/Tc  is  forced  to  equal  2$SpCUgp,  by  the  use  of  sufficiently  high 
gain,  the  system  saS  be  driven  unstable.  To  better  appreciate  this  possi¬ 
bility,  consider  the  root  locus  plot  given  in  Fig.  4.  The  relationship  of 
Eq  4  is  graphically  Illustrated  here,  and  conclusions  as  to  the  maximum 
value  of  l/Tc  and  minimum  value  of  2(^0^  are  clearly  shown  to  be 


Probable  values  for  the  right  6ide  of  Eq  5  can  be  obtained  by 
considering  the  approximate  airframe  transfer  function  factors  given  in 
Ref.  17, 

2?sp“8p  -  -(2v  +  Mq  +  Mi) 

1  .  ^5 

W2  - 
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(6) 


AC* "*  *  j 


I  Denotes  dosed -loop  roots  for 
given  value  of  forward  loop  gain. 


Kp 

Kr 

X 

d 

Forward  Loop  Gain  =  (Kp)(KF)(Ag) 


Figure  4.  Locus  of  Closed-Loop  Roots 
for  Pitch  Attitude  Control 


?.P. 


(7) 


1  .  ^6e 

so  that  2{sp%p  -  ji-  = 

.  pSU0c2  /  2k?  \ 

where  1§  =  cC^g/Ci^  is  the  effective  elevator  control  arm  (measured  posi- 
tive  forward) .  For  conventional  tail-aft  airplanes  with  some  static 
margin  (lg  and  both  negative)  the  bracketed  terms  of  Eq  7  are  always 
positive.  For  a  canard  control,  lg  will  change  sign,  and  the  contri¬ 
bution  will  be  negative,  so  it  is  conceivable  that  the  entire  term  could 

i 

be  negative.  However,  a  general  observation  is  that  it  will  take  a  very 
unusual  kind  of  configuration  (small  lg,  low  values  of  -Cm0>  which,  for 
conventional  airplanes,  is  usually  an  order  of  magnitude  greater  than 
C}^,  etc.)  to  make  the  value  of  25SpOjgp  -  l/T g2  negative,  as  sketched  in 
Fig.  4.  Therefore  only  for  such  unusual  airframe  configurations  is  there 
a  possibility  of  driving  to  zero  to  achieve  a  sustained  PIO. 

The  foregoing  demonstrates  that  for  airplanes  with  negligible  control 
system  dynamics  (including  nonlinear  elements  or  bobweight  effects), 
longitudinal  PIO's  involving  only  attitude  control  are  essentially 
impossible.  However,  in  variable- stability  flight  testing  and  ground- 
based  simulation  studies  where  the  general  practice  is  to  hold  Zy  constant 
and  vary  £Sp,  o^p,  stick  force  and  displacement  per  g,  etc.,  artificial 
relationships  between  £sp<10cp  and  l/T can  lead  to  PIO's  of  the  simple 
type  under  consideration  here.  For  example,  the  data  of  Ref.  5>  repro¬ 
duced  in  Fig.  5,  show  PIO  "tendencies"  for  the  lower  left  region  of  the 
£,  enplane.  These  data  were  obtained  for  a  fixed  value  of  Zy  correspond¬ 
ing  to  l/T@2  =  5*22  for  all  the  conditions  tested.  The  theoretical 
boundary  for  zero  £g0  as  given  in  Eq  is  superposed  on  the  original  plot. 
It  may  be  seen  tnat  there  is  general  agreement  between  the  predicted 
possible  "simple"  PIO  region  to  the  left  of  the  boundary  and  the  observed 
region.  The  fact  that  the  experimental  region  for  very  light  stick-force 

gradients  lies  somewhat  to  the  right  of  the  theoretical  boundary  is 

« 

evidence  of  additional  dynamics — in  this  case  nonlinear  effects  due  to 
the  high  ratio  of  breakout-force- to- stick-gradient,  1 .2  Ib/l .0  lb/g. 
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j  1  „  ;  Although  che  possibility  of  simple  PIO's  exists  to  the  left  of  the 

j  boundary,  whether  they  in  fact  will  occur  depends  on  the  likely  value  of 

l  pilot-vehicle  system  gain.  Thus  a  further  analysis  requires  an  examina- 

! 

!  tion  of  the  probable  pilot's  "mid-frequency"  gain  for  normal  closed- loop 

•  control  prior  to  the  PIO,  and  comparison  of  this  gain  level  with  the 

mid-frequency  gain  required  to  sustain  an  oscillation.  The  process  is 

j 

‘  illustrated  in  Fig.  6.  Here  the  condition  chosen  corresponds  to  the 

|  lower  left  data  point  in  Fig.  5  (f  *=  0.5  cps,  £  =  0.2).  The  neutral 

;  stability  point  marked  on  the  solid  bode  (corresponding  to  -180°  phase) 

J 

|  represents  the  pilot-airframe  combination  for  a  sustained  oscillation, 

j  The  dynamics  indicated  on  this  Bode  plot  are  all  due  to  the  airplane, 

since  the  pilot's  describing  function  is  simply  a  gain.  The  dotted  Bode 
represents  conditions  for  closed-loop  piloted  control  prior  to  the  PIO. 
Here  the  pilot  has  adapted  moderate  lead  in  accordance  with  the  closed- 
loop  adaptation  rules  of  Fig.  2.*  The  value  of  i/Tl  should  be  near  <%p, 
and  for  convenience  has  been  selected  equal  to  it.  Also,  as  always  for 
compensatory  situations  with  random-appearing  inputs,  the  pilot  model 
includes  the  reaction  time  delay  term,  e“Ts.  The  dotted  gain  line  shown 
as  appropriate  to  closed-loop  control  is  set  to  give  a  phase  margin  of 
about  40°  and  a  gain  margin  of  about  6  db.  Further,  the  asymptotic 
crossover  frequency  (essentially  the  closed-Jpop  bandwidth)  is  about 
1  rad/sec  consistent  with  usual  adjustment  criteria. 

The  important  point  in  looking  at  these  two  plots  is  the  gain  change 
required  to  go  from  a  conpensatory  tracking  to  an  oscillatory  (PIO)  con¬ 
trol  situation.  In  this  case  an  increase  of  about  a  factor  of  two  is 
required.  This  is  not  a  large  change  as  regards  pilot  adaptability 
provided  the  stick  forces  involved  are  not  excessive.  Thus,  for  suffi¬ 
ciently  low  stick-force  gradients  PIO  is  likely.  Conversely,  as  the  force 
gradient  increases,  the  physical  effort  required  to  increase  gain  by  a 
factor  of  two  serves  either  to  warn  the  pilot  against  this  course  of 
action  or  to  completely  prevent  it  (if  the  forces  ore  very  large). 


*Incidentially,  it  is  observed  in  Ref.  5  that  pilot  behavior  in  the 
"toe"  region  of  Fig.  5  resembled  a  lead. 
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Figure  6.  Bode  Diagram 
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.A'l%6ugfe/-^he.-^I0  region  of/ Mg.  §  is  thus  shown  to  be  dependent  on 
-a^tifltci^l  -re'l^tion'sldp'B-  tefcvijeen  1  |Te2;  tsp<%p>  there  att  flight  test 
examples  of  PIO's  in  precisely  the  same  region  (e.g.,  Refs.  6  and  18) . 


The  degree  to  which  all  such  situations  are  dependent  on  control  system 


or  bobweight  contributions,  linear  or  nonlinear,  is  not  precisely  known. 
Nevertheless  it  appears  to  be  true  that  longitudinal  PIO’s  were  non¬ 
existent  (or  not  reported)  until  the  advent  of  modern  hydraulically 
powered  elevator  actuation  systems.  The  PIO  of  Ref.  ’8  can  in  fact  be 
traced  directly  to  the  linear  contribution  of  the  hydraulic  system.  That 
is,  an  analysis  similar  to  that  of  Fig.  6  (including  the  linear  contribu¬ 
tion  of  the  hydraulic  system  as  measured  and  reported  in  Ref.  18)  shows 
that  the  complete  system  is  unstable  at  a  frequency  of  about  5*4  rad/sec. 
The  actual  PIO  record  shows  a  basic  frequency  near  6  rad/sec. 


As  a  final  observation  it  may  be  concluded  that,  except  for  very 
unusual  configurations,  longitudinal  attitude  control  PIO’s  can  be  sus¬ 
tained  only  for  conditions  in  which  control  system  dynamics  are  a 
contributing  cause. 


2.  Exaaple  2:  PIO  During  Roll  Control 


For  control  of  bank  angle  with  aileron,  the  pilot-vehicle  open-loop 
describing  function  for  potential  PIO's  is  given  as  follows,  remembering 
that  Yp  A  Kt»  during  a  pure  oscillation: 


Yp(s)  f  (r) 


KP^a(s2  +  2^S  +  4) 

(8  +  (6  +  ^r)(s2  +  2^B  + 


(8) 


t 

As  noted  in  earlier  theoretical  work  (Refs.  9,  12,  and  19)  there  can  be  a 
strong  effect"  tending  to  produce  an  oscillatory  instability  near 

Dutch  roll  frequency.  Recent  flight  test  examinations  (Ref.  20)  have 
confirmed  this.  The  effect  is  easily  seen  on  the  generic  root  locus  of 
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Fig.  7  which  shows  that  >  1  and  relatively  low  values  of  £<3,,  £<p> 

and  are  the  basic  requirements  for  a  sustained  oscillation.  Whether 
or  not  such  oscillations  occur  is,  again,  dependent  on  the  gain  differ¬ 
ences  between  compensatory  closed-loop  control  and  that  required  to  drive 
the  system  unstable.  Thus,  for  situations  in  which  the  o^,  pair  are 
well  separated  from  the  usual  crossover  frequency  region  (desired  closed- 
loop  bandwidth)  the  probability  of  a  PIO  is  low,  unless  the  daubing  ratio 
is  very  small  ( £  <  0.05)  as  for  some  high  altitude  aircraft.  A  lateral 
PIO  can  also  result  if  some  unusual  forcing  function  requires  a  change  in 
the  crossover  frequency.  For  a  more  complete  discussion  of  the  complex 
factors  involved,  see  Ref.  12,  13  and  20. 

Alternate  piloting  techniques  can  reduce  the  probability  of  these 
PIO's.  For  example,  use  of  the  rud*der  to  damp  the  Dutch  roll  will  increase 
both  and  The  effect  is  to  move  the  looped  root  locus  connecting 
with  ufy  well  into  the  (stable)  left  half -plane.  Or  the  pilot  may  learn 
to  fly  at  a  suppressed  gain  level  to  avoid  exciting  the  oscillation.  If 
such  techniques  (or  others)  are  required  to  cope  with  the  situation,  they 
will  not  be  generally  acceptable  to  pilots  and  will  be  construed  as  bad 
handling  qualities  (Ref.  21),  Furthermore,  since  they  are  not  natural 
and  instinctive,  they  cannot  be  relied  upon  to  prevent  PIO's  in  emergency 
conditions. 

The  basic  cure  for  such  PlO-prone  situations  involves  either  a 

>  reduction  in  ofy/ay  (by  cross-feeding  aileron  to  rudder,  or  by  suitable 

augmentation  driving  into  -  e.g.,  large  roll  damping)  or  an  increase 

% 

in  ^  and  (by  proper  feedbacks  to  a  yaw  damper) . 

3.  Example  3:  Rata-Baturatad  Elevator 

This  example  considers  a  PIO  resulting  when  the  pilot's  stick  movements 
are  faster  than  the  maximum  surface  rate  available  from  a  ‘surface  servo- 
actuator.  In  a  hydraulic  positional  servo,  subjected  to  periodic  inputs 
of  this  nature,  flow-rate-limiting  leads  to  a  reduction  in  the  ratio  of 
surface  velocity  amplitude  to  the  amplitude  of  the  servo's  oscillating 
,  error  signal,  thereby  reducing  the  forward  loop  gain  of  the  positional 

& 


Figure  7*  Root  Locus  Diagram  i'cr  the  uty/uy  Efl’ect 
During  Roll  Control 
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servo  and  thus  increasing  its  effective  tine  constant.  Qualitatively,  the 
result  is  a  narked  degradation  in  the  servo  characteristics,  effectively 
shifting  its  dynamics  to  lower  frequencies.  The  major  effect  is  an 
increase  in  the  phase  lag  at  frequencies  near  the  pilot-vehicle  system 
crossover.  This  effect,  while  somewhat  offset  by  an  associated  attenua¬ 
tion  increase,  can  still  sometimes  be  sufficient  to  result  in  a  limit 
cycle  when  large  surface  commands  are  present. 

The  analysis  of  such  series  nonlinearities  is  best  done  by  describing 
function  techniques,  previously  introduced  in  Section  II,-  but  in  a 
graphical  rather  than  analytical  form. 

The  basic  concept  of  describing  function  analysis  as  applied  to 
prediction  of  closed- loop  limit  cycles  (PIO)  is  to  assume  that  so6e  limit 
cycle  exists,  replace  the  actual  waveforms  around  the  loop  with  their 
Fourier  fundamentals,  and  then  to  see  if  conditions  are  satisfied  for  the 
oscillation  to  be  sustained.  The  criterion  for  a  neutrally  damped  oscilla¬ 
tion  is  simply  that  the  open-loop  amplitude  ratio  is  1.0  and  the  phase  —180°. 

When  the  nonlinear  elements  are  in  series,  the  system  dynamics  are 
separated  into  (a)  a  linear  portion,  represented  by  a  frequency-dependent 
transfer  function,  G(s)  [or  more  properly  G(jo>),  since  only  sinusoids  are 
considered] }  and  (b)  a  nonlinear  portion  which  has  transfer  characteristics 
represented,  in  general,  by  a  frequency-  and  amplitude-dependent  describing 
function j  N(j<n,  A/a).  (A/a  represents  the  signal  amplitude  relative  to  the 
nonlinearity.)  Then  for  an  oscillation  persist, 

G(jou)  •  N(ju),  A/e)  ■  -1 
or 

One  way  of  using  the  describing  function  technique  is  to  compute  or 
measure  the  frequency  response  of  the  nonlinear  element  for  a  set  of 
constant- input  amplitudes,  and  use  a  conventional  Bode  plot  to  plot  the 
total  open- loop  frequency  response,  Yq^  *  G  •  N,  for  each  member  of  the 
set.  Then  the  locus  of  intersections  of  the  unity-amplitude-ratio  points 


(9) 

do) 
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(0  db  crossovers)  is  projected  to  each  corresponding  phase  curve,  and  the 
locus  of  phase  margins  is  found.  If  the  phase  margin  locus  crosses  the 
180°  line,  then  a  limit  cycle  can  occur  at  the  frequency  and  level  of 
input  corresponding  to  the  intersection.  The  stability  of  this  limit 
•cycle  is  readily  determined  as  follows: 

The  limit  cycle  will  be  stable  (persist)  if  the  phase  or 
gain  margins  become  more  positive  for  increasing  oscilla¬ 
tion  amplitudes,  and  unstable  (die  out  or  diverge)  if  the 
phase  or  gain  margins  become  negative  for  increasing 
oscillation  amplitudes. 

A  particularly  single  describing  function  technique  is  possible  when 
the  describing  function  is  dependent  only  on  oscillation  amplitude,  and  is 
independent  of  frequency.  Making  use  of  Eq  10,  the  negative  inverse 
describing  function  -l/N(A/a)  is  plotted  on  a  standard  rectangular-grid 
gain  versus  phase  plot  with  A/a  as  a  parameter.  The  linear  portion,  G(ju>), 
is  also  plotted  on  the  gain-phase  plot,  but  in  the  conventional  fashion 
(i.e.,  using  frequency  as  a  parameter)  appropriate  to  Nichols'  chart 
analysis.  An  intersection  of  these  two  curves  satisfies  Eq  10  and  gives 
the  frequency  and  amplitude  of  the  limit  cycle.  The  stability  of  the 
limit  cycle  i6  given  by  the  same  criterion  as  before,  i.e.,  positive  gain 
or  phase  margins  for  increasing  oscillation  amplitude.  Thi6  technique  is 
used  for  the  case  at  hand  (see  Fig.  8). 

The  linear  transfer  function  plotted  in  Fig.  8  is  that  corresponding 
to  elevator  control  of  pitch  attitude  (Eq  1 )  for  a  PIO  encountered  on 
landing  of  the  X-15  airplane  (Ref.  22).  The  value  used  for  the  short- 
period  frequency,  <%p  i  2.5,  is  based  on  an  observed  "elevator-fixed" 
oscillation  occurring  at, about  the  same  IAS  as  that  for  PIO  (but  about 
80  sec  prior  to  PIO  onset).  The  values  of  l/T g2  and  ^sp^p  were  esti¬ 
mated  from  the  basic  aerodynamic  data  of  Ref.  25  as  0.82  and  1 .42, 
respectively.  The  gain  used  in  plotting  the  o/5e  transfer  function  of 
Fig.  8  corresponds  to  a  good  linear  loop  closure. 

Also  plotted  in  Fig.  8  is  the  inverse  of  a  simplified  approximation 
to  the  describing  function  for  the  rate- limited  positional  servo  (see 
the  Appendix  for  details  of  its  development).  Two  intersections  of  -1 /n 
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and  G  are  present  for  the  gain  used  in  plotting  the  linear  elements; 
a  stable  limit  cycle  is  indicated  at  to  «  2,  and  an  unstable  limit  cycle 
is  at  to  =  3‘  If  the  loop  gain  is  slightly  reduced  (i.e.,  the  entire 
linear  plot  shifted  down) ,  it  would  show  that  at  the  tangent  point  a  limit 
cycle  would  start  and  persist  at  about  2.8  rad/sec,  when  the  pilot  commanded 
an  elevator  rate,  6C,  about  1 .5  times  the  saturation  rate,  §s.  Since  the 
elevator  rate  is  saturated  for  most  of  each  cycle,  the  magnitude  of  the 
elevator  oscillation  is  given  by: 

l»l  -  if  O') 

For  the  limiting  elevator  rate  of  15°/sec  used  on  the  first  X-15  flight 
(Ref.  24)  and  the  theoretical  frequency  of  2.8  rad/sec  determined  above, 
the  expected  maximum  value  of  Be  would  be  ±5*3°  for  a  total  of  10.6°  peak 
to  peak.  The  flight  test  records  show  that  at  the  inception  of  PIO  (where 
the  aerodynamics  were  evaluated)  the  period  of  the  oscillation  was  about 
2  sec  (o>  =  3  rad/sec)  and  the  total  elevator  excursion  was  about  11°. 

This  agreement  between  predicted  and  actual  behavior  is  perhaps  fortuitous 
in  view  of  the  simplifying  assumptions  made  in  the  nonlinear  analysis  given 
in  the  Appendix.  Nevertheless  the  analysis  does  permit  a  basic  understand¬ 
ing  of  and  appreciation  for  the  design  implications  of  the  problem. 

For  example,  increasing  the  surface  velocity  limits  does  not  necessarily 
guarantee  that  PIO’s  will  be  eliminated.  A  change  in  6S  leaves  the  normal¬ 
ized  -1  /N  curves  of  Fig.  8  completely  unaltered  and  simply  implies  that  to 
get  PIO* 8  the  pilot  must  new  move  the  stick  at  a  correspondingly  increased 
maximum  rate.  Since  the  incipient  limit  cycle  frequency  is  also  unchanged, 
this  increased  rate  requires  a  larger  stick  position  input  and  results  in 
corresponding  increases  in  the  elevator  deflection  and  airplane  motions. 

If  this  rate  of  stick  motion  is  difficult  to  attain,  either  inherently  or 
because  of  the  high  stick  forces  required  (the  maximum  force  required  to 
move  the  stick  sinusoidally  against  a  bottoming  valve  depends  on  the 
flexibility  between  the  stick  and  valve),  then  the  increased  surface  rate 
will  probably  eliminate  PIO  tendencies.  On  the  other  hand,  if  the  increased 
stick  rates  are  easily  attained,  then  the  PIO  may  be  worse,  i.e.,  of  larger 
amplitude! 
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Increasing  the  short-period  damping  to  0.?>  and  again  keeping  the 
linear  open- loop  gain  con^atible  with  closed- loop  considerations,  gives 
the  dashed  gain-phase  plot  of  Fig.  8.  Not  until  the  pilot’s  gain  is 
raised  about  8  db  (a  factor  of  2.5)  will  these  modified  linear  character¬ 
istics  intersect  the  inverse  describing  function  plot  in  limit  cycle  con¬ 
ditions  described  by  a  frequency  of  about  2.6  rad/sec  and  an  5C/6S  =  2.0. 
Furthermore,  to  sustain  an  oscillation  requires  the  pilot  to  increase  the 
amplitude  of  his  input  motions  over  those  for  the  lower  damping  case  by 
the  ratios  of  5c/o,  a  factor  of  (2. 0/2.6) (2.8/l .5)  *  1.44.  Both  these 
influences  alleviate  the  PIO  susceptibility  of  the  system,  especially  if 
stick-force  gradients  are  of  reasonable  magnitude. 

As  a  matter  of  record,  the  fix  adopted  for  the  X-15  PIO  problem  took 
an  approach  combining  the  various  possibilities  outlined  above.  Quoting 
from  Ref.  24,  "...it  was  decided  by  the  manufacturer  that  (l)  the  less 
sensitive  center  stick  would  be  used  in  subsequent  landings,  (2)  the 
control- surface  rate  would  be  increased  from  15  deg/sec  to  25  deg/sec, 

(3)  the  longitudinal-force  gradient  would  be  increased  approximately 
30  percent,  (4)  the  longitudinal  breakout  force  would  be  increased 
slightly,  and,  as  an  additional  precaution,  (5)  launches  would  be  performed 
only  if  the  pitch  damper  were  operating." 

4.  Examp la  4:  Monllnaar  Bobveight-Frlction  Effects 

This  example  is  taken  directly  from  Ref.  1,  which  can  be  consulted  for 
more  detail.  It  is  included  in  the  present  report  as  an  illustration 
treating  a  fairly  complicated  feedback  nonlinearity  using  two  sets  of 
linearised  closed-loop  behavior,  which  correspond  to  the  two  extremes 
resulting  from  negligible  or  full  effect  of  a  subsidiary  loop  nonlinearity. 

Controlled  Element  Characteristic* .  For  an  airplane  equipped  with  a 
bobweight,  the  presence  of  friction  in  the  control  system  makes  the  air¬ 
plane’s  frequency  response  vary  as  a  function  of  input  force  amplitude. 

For  small  stick  forces  resulting  in  small  aircraft  accelerations,  the 
bobweight  feedback  fo**ce  is  too  small  to  break  through  the  friction  band; 
for  high  forces  and  resulting  accelerations  the  bobweight  ir  full  effective 


40 


while  for  intermediate  pilot  force  levels  the  bobweight  effect  varies 
between  inactive  and  active  states,  perhaps  in  rapid  succession.  Consider 
as  the  two  extremes  the  effective  controlled  element  characteristics  with 
no  bobweight  effect  and  with  full  bobweight  effect.  Kg  lb/g.  Then,  for 
the  system  depicted  in  the  block  diagrams  of  Figs.  3  and  9t  and  the  con¬ 
stants  given  in  Fig.  9,  the  resulting  0/Fs(s)  transfer  functions  are 
plotted  in  Bode  form  in  Figs.  10  and  11  (note  that  the  computation  of  the 
bobweight  case  requires  closing  the  inner  n£  loop).  The  high  altitude 
condition  (Fig.  10)  is  one  for  which  the  example  pilot-vehicle  combination 
never  encountered  PIO  problems;  for  the  low  altitude  condition  a  very 
severe  PIO,  of  about  ±8g  maximum  amplitude,  triggered  the  investigation  of 
Ref.  1 .  The  large  difference  in  the  controlled  element  dynamics  due  to 

vS 

the  variation  in  the  bobweight  influence  at  low  altitudes,  as  contrasted 
with  the  lesser  variation  at  high  altitudes,  is  +-he  primary  cause  of  the 
PIO  tendencies  in  this  case. 


Closed-Loop  Dynamic a.  The  analysis  behind  this  conclusion  starts  by 
considering  the  pilot-vehicle  combination.  For  the  low  altitude,  low  load 
factor  (zero  effective  bobweight)  condition,  the  pilot  adaptation  in  the 
frequency  region  of  interest  will  be  a  lag.  This  allows  him  to  get  good 
crossover  in  a  fairly  extensive  region  of  6  db/octave  slope  and  is  con¬ 
sistent  with  the  experimental  evidence  presented  in  Ref.  10.  In  the  case 
at  hand,  a  lag  time  constant  fairly  close  to  T it  desirable  and,  for 
convenience,  it  is  assumed  that  Tj  =  Tgg.  The  pilot  model  of  interest 
here  is  thus  simply 


-TS 


=  0.2 


(Txs  +  1 ) 


Ti 


«  tT  =  5‘18 

T02 


and  the  corresponding  phase  and  amplitude  contributions  result  in  the 
3olid  Bode  plot  of  Fig.  12.  Assuming  a  phase  margin  of  kO°  givoc  the 
probable  gain  crossover  shown,  resulting  also  in  a  gain  margin  of 
Km  =  4  db.  On  the  other  hand,  to  sustain  a  large  amplitude  PIO  (where 
the  bobweight  is  fully  active)  requires  a  pure  gain  level  corresponding 
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Figure  9*  Block  Diagram  and  Dynamics  for  Analysis 
of  Nonlinear  Bobveight  Friction 
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FLIGHT  CONDITIONS!  Ms 0.85  h  =  25000  ft 


Figure  10.  Effect  of  Bobweight  on  Pitch  Response 
Away  from  PIO  Conditions 


Amplitude  Ratio  (You|  * 


to  the  zero  phase  margin  of  the  dashed  Bode.,  This  gain,  level  is  between; 

7  and  12  db  less  than  that  for  compensatory  controli  Notice  also  that 
higher  gains — that  is ,  gains  closer  to  the  compensatory  level— will  produce 
diverging  oscillations.  It  appears,  therefore,  that  to  avoid  large  ampli- 
^  tude  PIO’s  the  pilot  must  reduce  his  gain  at  lease  by  a  factor  of  about 

3  to  4  from  that  normally  desired  (and  probably  used)  for  small  precise 
corrections.  The  excessively  nonlinear  and  rapid  readaptation  required  of 
the  pilot  to  avoid  PIO  is  more  than  the  pilot  can  accommodate.  Also,  the 
basic  characteristics  with  the  bobweight  active  (o^p  =  9.8,  £gp  =  0.10, 
l/Tg2  =  3*1 8)  are  themselves  quite  poor  (Ref.  5),  and  even  if  there  were 
no  nonlinear  action  due  to  friction  the  system  would  probably  be  PIO  prone 
in  the  sense  of  Exangle  1  (i.e.,  S^gpO^p  -  1/Tq2  =  -1.2). 

Pure.  The  basic  cure  lies  in  removing  the  peak  in  the  bobweight- 
augmented  response  while  preserving  the  low  frequency  3tick-force-per-g 
,  characteristics  desired  at  this  flight  condition.  This  is  done  most 

simply  by  reducing  the  bobweight  loop  gain  by  a  factor  of  2  to  4  while 
preserving  low  frequency  stick-force-per-g  characteristics  required  by 
specifications.  This  is  feasible  because  the  pilot  loop  gain  involves 
the  quotient  of  Kg/Kp  (where  Kp  is  the  feel  spring  gradient:  Kp  «  Fs/Sj^), 
while  the  static  stick  force  per  g  involves  the  sum  of  Kg  +  (Kp  *  constant) . 
The  bobweight  gain  was  reduced  from  Kg  =  2.0  to  1 .0  lb/g,  while  the  feel 
spring  gradient  was  increased  (for  small  stick  deflections)  from 
Kp  =  5  to  10  lb/in.  The  resulting  loop  gains  and  stick-force-per-g 
parameters  are  shown  in  Table  III.  WhereaB  the  original  T-38A  bobweight 


TABLE  III 

COMPARISON  OF  BASIC  AND  MODIFIED  PARAMETERS 

Basic  T-38A  Modified  T-38A 


Bobweight  gain,  lb/g . 

0 

2.0  (max) 

1 .0  (max) 

Feel  spring  gradient,  lb/in . 

5.0 

5-0 

10.0 

Gearing,  deg/in . 

1 .03 

1 .03 

1 .05 

Inner  loop  gain  at  PIO  conditions, 

Vi's . 

0 

1 .0 

0.25 

Static  stick  force  per  g,  lb/g . 

2.0 

4.0 

4.6 

Minimum  dynamic  stick  force  per  g, 
lb/g . 

1  .2 

0.7 

2.9 

4b 


characteristics  resulted  in  a  drop-off  in  stick  force  per  g  from  k.O  to 

0.7  lb/g  near  the  PIO  frequency,  the  modified  characteristics  only  drop 
from  4.6  to  2.9  lb/g  in  this  frequency  region.  It  is  also  apparent  that 
the  original  unaugmented  airframe  had  overly  sensitive  static  stick-force- 
per-g  characteristics.  It  is  felt  that  the  bobweight  loop  gain  redaction 
of  400  percent  is  about  the  minimum  which  will  result  in  elimination  of 
the  basic  PIO  tendencies,  although  lesser  percentages  would  certainly 
result  in  substantial  improvement . 

The  consequences  of  increasing  the  feel  spring  gradient  in  conjunction 
with  reducing  the  bobweight  gain  yield  a  twofold  benefit: 

1 .  The  bobweight  loop  gain  is  reduced  as  Kp  is 
increased. 

2.  The  feel  system  frequency,  oujr,  is  increased  by  the 
square  root  of  the  change  in  Kp,  i.e.,  from  18  to 
25  rad/sec,  which  moves  it  farther  away  from  the 
short-period  frequency  and  reduces  its  adverse 
effect  on  the  coupled  short-period  damping  ratio. 

Increasing  Kp  has  some  serious  disadvantages ,  however.  Since  the  bobweight 
influence  is  reduced,  the  static  stick  force  per  g  is  more  affected  by 
variations  in  the  flight  conditions,  so  Fs/g  will  not  be  as  constant  over 
the  entire  flight  envelope  as  in  the  basic  design.  Furthermore,  tne 
increased  force  gradient  may  result  in  excessive  maximum  stick  forces 
during  landing,  takeoff,  or  supersonic  conditions  when  large  stick  deflec¬ 
tions  are  required.  Thus,  the  cure  for  the  PIO  tendencies  will  perhaps 
detract  from  the  general  airplane  handling  qualities  taken  across  the 
entire  flight  regime. 

Optimum  Recovery  Technique.  This  analysis  also  reveals  the  optimum 
technique  for  rapid  recovery  from  a  pilot-induced  oscillation  of  this 
type.  If  a  PIO  is  started,  it  can  be  stopped  either  by  releasing  the 
stick  or  by  clamping  it  securely,  i.e.,  opening  the  loop.  Clamping  the 
stick  will  cause  the  oscillation  to  decay  at  the  stick-fixed  short-period 
damping  ratio  (£Sp  =  0.4),  whereas  releasing  the  stick  will  cause  the 
oscillation  to  decay  at  the  stick-free  short-period  damping  ratio  (£gp  =  0.1). 
Since  is  always  less  than  (;sp>  the  best  recovery  procedure  theoretically 
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1b  to  clamp  the  stick  securely*  For  the  unmodified.  T-j8A>  this  should 
damp  the  oscillations  roughly  2  to  4  times  as  fast  as  releasing  the  stick, 
depending  on  the  amount  of  stick  friction.  However,  if  it  is  difficult 
to  achieve  rigid  clanging  during  a  violent  PIO,  then  releasing  the  stick 
is  the  only  alternative. 

Summary.  In  summary,  it  is  concluded  that  at  low  altitude,  high 
subsonic  flight  conditions  the  airframe  sensitivity  to  stick  forces  is 
sufficiently  high  that  bobveight  effects  can  produce  marginally  stable 
stick- free  short-period  characteristics  at  large  load  factors.  At  low 
load  factors  control  system  friction  prevents  bobweight  feedback.  When 
large  stick  inputs  result  in  high  load  factors,  the  full  effect  of  the 
bobweight  interactions  are  felt,  and  the  pilot  must  reduce. his  gain 
severalfold  to  avoid  a  pilot-induced  oscillation.  To  stop  a  developed 
PIO,  clanging  the  stick  is  theoretically  more  effective  than  simply  1 
releasing  the  stick.  At  high  altitude  conditions,  the  difference  in 
tolerable  pilot  gain  for  instability  due  to  bobweight  effects  is  negligible, 
and  no  PIO  is  induced.  The  simplest  recommended  cure  consists  of  increas¬ 
ing  the  feel  system  spring  gradient  for  small  stick  deflections  while 
decreasing  the  bobweight  gain. 

j}.  Oonaludlng  Remarks 

The  wide  range  of  PIO  problem  areas,  physical  causes,  analysis  techniques, 
and  cures  illustrated  in  the  preceding  examples  shows  that  pat  prescriptions 
for  PIO  causes,  analyses,  and  cures  are  not  to  be  expected.  Each  case  must 
be  treated  individually  and  thoroughly,  and  competing  fixes  should  be  evolved 
before  a  decision  as  to  the  best  is  made.  The  hopeful  note  in  all  this  com¬ 
plexity  is  the  demonstration  herein  that  well-known  analytical  techniques, 
combined  with  recent  data  on  the  human  pilot's  behavior  (both  prior  to  and 
during  a  PIO)  and  system  describing  function  data  in  the  proper  detail,  can 
correctly  assess  the  basic  causes  of  pilot-induced  oscillations  and  reveal 
the  most  promising  cures. 
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APPBDH 


S0XVASXCB  07  TEE  D»CRIB1JG  fUBCTIOB 
TOR  A  RATX-8A7URATID  P06ZTIORAL  SERVO 


The  rate-limited  servo  is  described  by  the  following  block  diagram 
and  parameters: 

Saturation 


The  output  rate  is  proportional  to  the  error  with  gain  K  for  errors 
less  than  €s,  but  is  saturated  at  a  level  &s  for  error  signals  larger 
than  es.  The  problem  in  determining  the  closed-loop  sinusoidal  describing 
function,  6/6c,  is  that  the  error  signal  is  not  sinusoidal  when  the  rate 
exceeds  saturation,  and  it  is  not  simply  related  to  the  input,  as  it  is 
for  linear  operation.  A  numerical  or  graphical  cross-plotting  scheme, 
such  as  that  described  in  Ref.  1  5  ,  pp.  209  ff»,  could  be  used  to  obtain  an 
"exact"  describing  function  for  &/5c.  A  simpler  analytical  approximation 
is  used  herein  which  reveals  all  the  essential  features  with- sufficient 
accuracy  for  most  practical  problems. 

The  error  signal  is  assumed  to  be  closely  represented  by  a  pure 
sinusoid  of  amplitude  E:  e  «  E  sin  (cut  +  cp).  The  sinusoidal  describing 
function  for  the  limiter  itself  is  well  known  (e.g.,  Ref.  19,  pp.  106, 

110,  11 4,  129,  258).  For  E  >  es  the  amplitudes  of  the  Fourier  fundamental 
of  the  output  is  given  in  terms  of  the  input  to  the  nonlinearity  by 


The  effective  gain  of  the 

nonlinearity  is 

• 

/ 

K‘ 

=  K 

for  E  £  s6 

(A-2) 

K* 

r? 

for  E  >  €s 

(A-5) 

This  relationship,  as  a  function  of  E/es,  is  shown  in  Fig.  A-1 ,  normalized 
with  respect  to  the  linear  gain,  K.  Also  shown  are  two  approximations  to 
the  saturated  describing  function: 


1. 


The  upper  approximation  is  asymptotically  correct  at 
large  error  amplitudes  where  the  output  rate  looks 
like  a  square  wave  (whose  fundamental  amplitude  is 
4/*  times  the  saturation  rate).  For  this  case, 


K* 


iia  .  E.»! 

«  E  *  €g  n 


(A-4) 


2.  The  lower  approximation  curve  is  asymptotically 
correct  at  errors  just  exceeding  the  saturation 
point,  where  the  slightly  clipped  output  still  looks 
like  a  sinusoid  whose  amplitude  is  approximated  by 
the  saturation  rate  itself,  6S,  instead  of  the 
Fourier  fundamental  of  Eq  A-1 . 
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6S 


1  .0 
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(A-5) 


It  is  apparent  from  Fig.  A-1  that  the  approximation  of  Eq  A-4  is 
actually  very  good  for  E/es  >  2.0,  and  is  about  as  good  as  Eq  A-5  for 
1 .0+  <  E/ €s  <  2.0.  Nevertheless,  Eq  A-5  is  simpler  to  handle  and  still 
reveals  the  essential  features  of  the  closed-loop  system. 


The  error  amplitude,  E,  must  now  be  related  to  the  command,  6C.  From 
fundamental  closed-loop  relationships  we  have 


Ye(s) 


6(b)  _  1  .  Y  _  K1 

§i ( s )  1  +  Yol(s)  *  0L  -  s 


(A-*) 


and  replacing  s  by  j<o  for  a  pure  sinusoidal  input  gives 


Carrying  out  the  algebra,  remembering  that  «  -i ,  and  solving  for  the 
magnitude  of  the  complex  quantities,  ve  get 


je(ja))l  J2_  co 

™  '  6= '  WvfF 


(A-7) 


where  6C  is  the  maximum  command  displacement  input*  At  this  point  K*  is 
still  dependent  on  E  in  general,  but  choosing  one  of  the  approximations 
of  Eq  A-4  or  A-5  removes  this  dependence  (for  E  >  e6)  and  is  the  analytical 
key  to  a  simple  solution.  Equation  A-7  may  then  be  manipulated  algebraically 
to  yield  expressions  for  |E/es|  and  thus  for  |e/Sb|  and  hence  K* ,  in  terms 
of  the  command  instead  of  the  error.  The  results  of  this,  using  the  simpler 
expressions  of  Eq  A-5,  are: 


The  closed-loop  response  can  now  be  determined  from  the  basic  relation¬ 
ship  for  YqL: 
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y0L  1 

1  +  Yot  “  1 
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Ts  +  1 


(A-12) 


For  the  linear  range  of  operation  (E/es  <  1.0)  this  is  just  a  first-order 
lag  given  as  T  =  1 /K,  which  increases  as  K1  becomes  less  than  K  beyond 
the  saturation  level. 
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Putting  s  =  jco  and  substituting  Eq  A-11  for  1  /K*  -gives 


Thus,  at  constant  input  amplitude,  the  closed-loop  describing  function 
depends  on  both  amplitude  and  frequency  once  saturation  occurs.  However, 
recognizing  that  the  input  rate  for  a  sinusoid  is  just  Sc  =  6cco,  Eq  A-l 5 
may  be  simplified  further  to  eliminate  the  frequency  dependence: 


(A-i  b) 


The  associated  gain  and  phase  are  given  by 


Since  Eq  A-9  is  asymptotically  correct  at  the  saturation  point,  we  can 
solve  it  for  %  at  saturation  (E/es  *=  1.0),  giving 


<«% 
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(A-l  7) 


which,  by  further  algebra  and  with  Sc  =  6co),  may  be  put  in  the  form 


2 

s 


>  1 


(A-18) 


The  meaning  of  this  is  as  follows:  K  is  the  linear  closed- loop  break 
frequency  of  the  servo  and  is  a  large  value  like  20  to  50  rad/sec.  The 
input  frequencies,  <n,  are  rather  low,  usually  less  than  10  rad/sec. 

Hence,  o^/K^  «  1  for  practical  problems,  and  the  saturation  criterion 
reduces  to  the  simple  statement  that  the  maximum  commanded  rate  exceeds 
the  output  rate  limit.  This  is  not  surprising,  since  the  output  is 
closely  following  the  commands  at  low  frequencies;  hence  the  output  rate 
starts  saturating  when  the  command  rate  reaches  the  saturation  value. 
Equation  A-18  shows  that  at  higher  frequencies  (comparable  with  the  servo 
break  frequency,  K)  more  command  rate  is  required  to  saturate  the  servo. 

Expressions  similar  to  Eq  A-8  through  A-1 8  result  for  the  asymptotic 
approximation  of  Eq  A-4,  with  the  constant  k/n  -  1.27  appearing  in  several 
places.  All  possible  forms  are  summarized  below. 

Linear  range: 


*  yCL  =  "tan_1  Y 

Saturation :  • 

Satu^tion  frequency: 


(A-1 9) 


(A-20) 
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For  E/es  1 .0+  (6C/Ss  1 .0+) ;  Small  saturation: 


l%l  1  o  «  r 


1  -* 


-i  .//M 


r)  -  ’ 


.-J«? 


-tan  4/  To  ~  1 


(A-21) 


For  E/cs  »  1.0  (&C/6S  »  1.0);  Large  saturation: 


IyclI  “ 


«  Bc to  °r  n  £c 
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The  resulting  approximate  describing  functions,  N(jo),  Bc/66),  are 
given  in  Table  A-I  for  various  values  of  5C/6S,  while  in  Fig.  A-2  the 
negative  inverse  describing  functions,  — 1  /N,  are  plotted  on  a  log-gain- 
phase  plot.  The  gain  phase  diagram  (Nichols  Chart)  facilitates  limit 
cycle  analyses  when  the  describing  function  depends  only  on  a  relative 
amplitude  parameter,  and  not  on  frequency.  Actually,  the  complete 
describing  function  for  the  rate-limited  servo  consists  of  three  parts: 
a  linear  portion  below  saturation  which  depends  only  on  frequency 
(first-order  lag  with  break  frequency  »  K) ;  a  second  small  portion  which 
depends  on  amplitude,  &c/&s,  and  applies  for  E/es  =  1 .0+;  this  blends  into 
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'  TABLE  A-I 

DESCRIBING  FUNCTION  FOR  A  RATE-SATURATED  POSITION  SERVO 
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a  third  portion,  beyond  E/es  =  2,  also  dependent  only  on  amplitude.  As 
shown  in  Fig.  A-2,  when  plotted  on  a  gain  phase  diagram,  all  three  portions 
lie  along  one  curve.  (The  linear  portion  has  frequency  as  a  parameter 
with  the  break  frequency  <jq  =  l/T  *=  K  at  |— i/n[  =  3  db  and  $  — l/N  =  135°*) 

Neutral  stability  of  the  linear  portion  occurs  when  the  G(jci5)  curve 
intersects  the  -l/N  curve  at  tha  laaa  fraquancy,  and  it  will  be  assumed 
that  positive  gain  and  phase  margins  exist  for  linear  operation  (small 
commands).  The  nonlinear  portion  then  becomes,  to  first-order  accuracy, 
dependent  only  on  the  relative  amplitude  of  command  rate  to  saturation 


rate  of  the  position  se.n/o,  5C/6S.  As  noted  in  Pig-  A-2,  the  two  different 
approximations  shift  only  the  5c/bs  parameter  along  the  curve  a  small 
amount,  probably  within  the  experimental  accuracy  of  the  known  saturation 
rate  and  form  of  the  saturation  curve.  Accordingly,  it  is  recommended 
that  the  expressions  for  the  simpler  form  of  approximate  gain  given  in 
Eq  A-5  and  A-21  be  used  for  analytical  purposes. 
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